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1. Introduction 

The study of the low-lying eigenmodes of the Dirac operator (LDE), those corresponding to 
the smallest eigenvalues, has a rich history since these modes are thought to be representative of, if 
not responsible for, much of the infrared behavior in QCD. Such effects include 

• Chiral Symmetry Breaking a la Banks and Casher where (yy) ~ p x (0) 

• The low eigenmodes of Tp dominate quark propagators 

• Confinement in many scenarios is thought to be related to topological excitations: instantons, 
monopoles, or vortices. These objects all localize Dirac zero-modes in some way. 

Thus we are interested to learn what we can about this localization, if indeed it exists, and then 
to characterize it in some quantitative way. Hopefully this characterization can tell us something 
about the mechanisms responsible for localization. 

2. Inverse Participation Ratio 

The Inverse Participation Ratio (IPR) provides a quantitative number which characterizes the 
localization of a scalar field. For the LDEs, it is defined as 

4 = VJ> 2 (*) 

X 

where V is the number of lattice sites x, and i labels which eigenmode is under consideration, 
V/(jc) is the i-th lowest eigenvector of the (asqtad) Dirac operator and 

X 

The IPR takes the following values in cases of different localization: 
1=1 if p is constant, 

/ = 1// if p is localized (and constant) on a fraction / of sites, and 
I = V ifp = 5„ . 

As was first pointed out in [j|], the fraction of points involved in localization should scale with 
the dimension of the localizing manifold. For example, 

Sf/a . £f/a 2 

/l -dim = TJ-J-7 or /2-dim = 7} , 7 , etc. 
V far V/a^ 

where in the first, one-dimensional, case «Sf is the total length of "localizing material" and lattice 
spacing a, or in the two-dimensional case, srf is the area of two-dimensional material, and so on 
for higher dimensions. The (possibly fractal) dimension, d, of the localizing manifold is given by 
the scaling of the IPR as the lattice spacing is varied, 

I = l/f~a d - A (2.1) 
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In recent years a growing number of authors have found evidence for localization of either the 
LDE or other quantities, such as the topological charge density Since the appearance of [[lj] 
other groups have analyzed the scaling of the IPR using improved ^ and alternative [Q] operators, 
finding similar conclusions. 

While the scaling of the IPR is a rather clean indication of the underlying localization dimen- 
sion, its extraction requires a wide range of lattice spacings to get a fit with good statistics. In [p 
for example, the dimension was given as somewhere between 2 and 3. For best results, scaling 
measurements should be done at fixed physical volume, which can then be compared to similar 
measurements at fixed lattice spacing to understand finite size issues. 

3. Scaling of the IPR 

In [[I]] we presented preliminary results for the scaling of the IPR using quenched lattices 
(Symanzik 1-loop improved gauge action) ranging from 12 4 and a = 0.2 fm up to 24 4 and a ~ 0.1 
fm. However, with better statistics on the finest lattices, we found that the lattice spacing was 
closer to a = 0.095 fm, 5% from of our target of a = 0.1 fm. In order to maintain a fixed volume 
we regenerated this ensemble at the target lattice spacing. Furthermore, as our lattice spacing is 
set using r\ from the static quark potential, we use here an updated value of r\ (0.317 fm) taken 
from [Q] (n = 0.344 fm was used in [[T|]). This increases all lattice spacings by ~10%, but has 
no effect on our results. Additionally, we have increased the convergence criteria for computing 
eigenvectors and added a finer 28 4 ensemble. 

The lattices and parameters used for the present work are collected in the table below. 



a 


L 


vol 


/3 


no. configs. 


0.218 fm 


12 


(2.61) 4 fm 4 


7.56 


100 


0.163 


16 


(2.61) 4 


7.847 


100 


0.128 


20 


(2.56 f 


8.109 


100 


0.110 


24 


(2. 64) 4 


8 .295 


100 


. 0915 


28 


(2.56) 4 


8 . 527 


100 



with 64 eigenvectors per lattice. The italicized ensemble was rerun for better convergence of eigen- 
vectors (with very little change), while ensembles in teletype are new. 

Our main result is summarized in Figure 1, where we plot the scaling of the average IPR 
using the lowest 8 eigenvectors. Compared with we have better statistics on the IPR values 
and, particularly, the new value at a smaller lattice spacing (0.0915 fm) gives a much more precise 
scaling dimension for the IPR. 

The best fit to the data has d — 4 = 0.934 ± 0. 149 (see eq. 2. 1), thus the scaling dimension, d, 
of the localization manifold is essentially 3. 

4. Mobility Edge 

In analogy with the mechanism of Anderson localization in condensed matter physics, we can 
investigate whether our data shows a mobility edge, an energy above which quarks are delocalized 
and below which they show localization. This feature has been investigated by Golterman and 
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0.14 0.16 0.1 6 
Lattice spacing (Fermi) 



Figure 1: The scaling of the IPR for the lowest 8 eigenvectors, in red. In green is the best fit, and associated 
values. In blue is a fit to constant + constant/a for comparison. 

Shamir ^ and observed in SU(2) []|] using an overlap Dirac operator with exact zero modes. The 
signal is a reduction in the IPR (less localization) at some critical value for the eigenvalue-the 
mobility edge. 

In our SU(3) data the IPR values are rather low as compared to IPR values in SU(2) studies 
[@] (3-4 vs 5-20). One possible reason for this might be that the localization is due to topology of 
one of the SU(2) subgroups, while the other two subgroups randomize the eigenvectors. Whatever 
the reason, this is an interesting clue to understanding the localization. 

Only our 28 4 lattice at a = 0.0915 fm shows a weak indication of the mobility edge, shown in 
Figure 2. Its value in physical units, around 50 ~ 100 MeV is consistent with that seen in ^ . A 
study of the IPR on various volumes would be required to confirm this. 

5. Two-point Correlations 

While the IPR is a good quantitative indicator of localization, it only tells us the fraction of 
lattice sites where the eigenvector is large. If we rearranged the lattice sites we would obtain the 
same IPR. The two-point correlator on the other hand gives information on the connectedness of 
the eigenvector, and should drop off as ~ r^ -4 if the eigenvector is localized uniformly on a d- 
dimensional manifold. 

We have computed the "all-to-all" correlator, (p(x)p(y)) and present the average at different 
spatial separation \x — y\ for eigenvalue in Figure 3. These correlators have parity sawtooth 
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Figure 2: (IPR) vs eigenvalue in physical units. The mobility edge is weakly indicated by a decrease in IPR 
at about 60-100 MeV. 



behavior familiar in staggered fermion propagators, which lessens as the lattice spacing decreases. 
Taking the slope of the correlator from the even points, we show a line with d = 3.5 (a very similar 
result if achieved for the odd points which lie slightly below the line). We see that in physical units, 
the correlation disappears by about 1 fermi, for all of our lattices. 

6. Conclusions 

We have extended our study of the IPR of the lowest eigenvectors of the asqtad Dirac operator 
computed on quenched background gauge fields. A number of conclusions are apparent. 

• The IPR scaling ~ I /a implies a dimension d = 3 for the localizing manifold. The signal is 
quite a bit clearer than in our previous study. 

• We see a weak mobility edge in our finest lattices (only) at 50~100 MeV. This is consistent 
with other work. The fact that our signal is weak is likely attributed to the lack of exact zero 
modes of our Dirac operator and possibly to the larger gauge group (SU(3) vs SU(2)). 

• The two point correlator (p(x)p(x + r)) suggests a fractal dimension of ~ 3.5. It should be 
noted that we have not finished computing this quantity on all of our largest lattices (Figure 
3 represents only about 8 of our 28 4 dataset, however the correlator shows little fluctuation 
from lattice to lattice). In physical units the correlation disappears at ~1 fm. 

Our findings do not support the naive picture where the low-lying Dirac eigenmodes are lo- 
calized on monopoles (<i=l) or vortices (<i=2). The relationship between topological excitations 
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Figure 3: The two-point correlator of the lowest eigenvector, J^, (p (x)p(x + r)) as a function of distance r. 
Plots (except a=0.0915) are displaced upward for clarity. The axes are log-log. 

and eigenvector localization is more subtle. For instance, a set of center vortices gives eigenvectors 
localized weakly on the vortices themselves, and more strongly on their intersections [j/J]. 

References 

[1] C. Aubin et ai, The scaling dimension of low lying dirac eigenmodes and of the topological charge 
density, Nucl. Phys. (Prvc.Suppl.) 140 (626) [hep-lat/0410024] 

[2] I. Horvath et al. On the local structure of topological charge fluctuations in QCD, Phys. Rev. D67 
(011501) [hep-lat/0203027] 

[3] F.V. Gubarev et al., Evidence for fine tuning of fermionic modes in lattice gluodynamics 
[hep-lat/0505016] 

[4] J. Greensite et al., Localized eigenmodes of covariant Laplacians in the Yang-Mills vacuum, Phys. 
Rev. D71 (114507) [hep-lat/0504008] 

[5] C. Aubin et al., Light pseudoscalar decay constants, quark masses, and low energy constants from 
three-flavor lattice QCD, Phys. Rev. D70 (114501) [hep-lat/0407028] 

[6] M. Golterman and Y. Shamir, Localization in lattice QCD, Phys. Rev. D68 (074501) 

[hep-lat/0306002] 

[7] H. Reinhardt et al, Quark zero modes in intersecting vortex gauge fields, Phys Rev. 66 (085004) 

[hep-th/0203012] 



299/6 



